Itinerant magnetic excitations in antiferromagnetic CaFe2As2 
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Neutron scattering measurements of the magnetic excitations in single crystals of antiferromag- 
netic CaFe2As2 reveal steeply dispersive and well-defined spin waves up to an energy of ~ 100 
meV. Magnetic excitations above 100 meV and up to the maximum energy of 200 meV are however 
broader in energy and momentum than the experimental resolution. While the low energy modes 
can be fit to a Heisenberg model, the total spectrum cannot be described as arising from excitations 
of a local moment system. Ab-initio calculations of the dynamic magnetic susceptibility suggest that 
the high energy behavior is dominated by the damping of spin waves by particle-hole excitations. 



There are indications [H, H, 0] that superconductiv- 
ity (SC) in the iron arsenides (Fe-As) family may be 
driven by a magnetic pairing mechanism, the nature of 
which remains poorly understood. Recent inelastic neu- 
tron scattering experiments have uncovered a magnetic 
resonance feature in superconducting Bao.6Ko.4Fe2As2 
4|, Ba(Feo.92Coo.o8)2As2 [5] and Ba(Feo.95Nio.o5)2As2 
6|, that indeed suggests a close relationship between 
magnetism and superconductivity. The interpretation 
of these data as excitonic or spin wave-like [8] de- 
pends crucially on the nature of the Fe-As magnetism, 
in particular the dimensionality and degree of itinerancy, 
which we here explore through magnetic neutron scatter- 
ing from the parent compound CaFe2As2. 

Recent experiments probing long wavelength spin 
waves in antiferromagnetic Fe-As single crystals revealed 
very large spin wave velocities and anisotropic three- 
dimensional magnetism Another measurement 
on polycrystalline BaFe2As2 showed excitations as high 
as 100 meV [nj, but the powder averaging inevitably re- 
stricts the information that can be extracted from the 
data. In this Letter, we report neutron scattering mea- 
surements from single crystals of CaFe2As2 probing mag- 
netic excitations up to 200 meV. Steeply dispersive spin 
waves, consistent with previous measurements, were ob- 
served up to ~ 100 meV. Beyond 100 meV, and up to the 
maximum observed magnetic excitation of ~ 200 meV, 
damping is substantial and collective excitations are not 
clearly visible in the data. Linear response calculations 
of the dynamic susceptibility reveal that the high energy 
antiferromagnons are strongly affected by particle-hole 
excitations (Landau damping), indicating that the spin 
excitations have an itinerant nature over a large frac- 
tion of the Brillouin zone. A Heisenberg model can ac- 
count for magnetic excitations below 100 meV and was 
used to extract effective low energy exchange parame- 
ters and an energy dependent spin wave relaxation rate. 



The local moment (Heisenberg) description of magnetism 
in this class of materials is however, limited to small 
wave- vectors and low energies with substantial effects of 
electron-itinerancy beyond 100 meV. 

CaFe2As2 is a non-superconducting parent compound 
of the new family of iron arsenide based superconducting 
materials [l2[. CaFe2As2 undergoes a paramagnetic to 
antiferromagnetic (AFM) transition simultaneously with 
a tetragonal (I4/mmm) to orthorhombic (Fmmm) struc- 
tural transition below Tat = = 172 K [i2|. The or- 
thorhombic lattice parameters are a = 5.51 A, 6 = 5.45 
A, and c = 11.67 A at T = 10 K. The magnetic struc- 
ture is a collinear antiferromagnet with propagation vec- 
tor Q,AFM — (101)- Throughout this paper, wavevectors 
are specified in the orthorhombic reciprocal lattice. The 
Fe ordered moment is 0.80(5)/iB/Fe directed along the 
orthorhombic a axis [l^. 

The measurements were performed on the MAPS time- 
of-flight chopper spectrometer at the ISIS neutron facil- 
ity. The sample consisted of approximately 400 crystals 
grown from tin flux, with a total mass ~ 2 grams. These 
were co-aligned with an angular mosaic width of 1.5 de- 
grees and a [H{)L] scattering plane. Crystals were etched 
in concentrated hydrochloric acid to remove Sn flux from 
the surface. Magnetic Brillouin zone (BZ) centers are lo- 
cated at Qafm = (i^ = 2m + 1, = 2n, L = 2j + 1), 
m, n, j being integers. The sample was mounted on the 
cold finger of a closed-cycle refrigerator and cooled to 
T = 10.0(5) K for measurements at four incident ener- 
gies; Ei = 100, 180, 300, and 450 meV with the incident 
beam along [100] and [001]. 

The quantity of interest is the dynamical structure fac- 
tor 5'(Q,a;), which is related to the double differential 
cross-section via, 5'(Q,a;) = {ki/kf)d?cF / dVtfdEf . Here 
subscripts i and / label the initial and final neutron wave 
vector /c, 1] is the solid angle of scattering, and is the 
scattered neutron energy. Incoherent nuclear scattering 
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FIG. 1: Observed (left) and calculated (right) magnetic excita- 
tions in CaFe2As2 single crystals at temperature T = 10 K with an 
incident energy Ei = 450 meV. The projections are in the scatter- 
ing plane formed by the energy transfer axis and (if, 0, 0) direction 
and averaged over —0.1 < K < 0.1. An estimated background 
away from the magnetic zone center averaged over 1.8 < H < 2.0 
and —0.25 < K < 0.25 has been subtracted from the data. The 
intensities have been multiplied by the energy transfer for viewing 
purposes. Due to the fixed crystal orientation with incident beam 
along L, the L component of the wavevector varies with energy 
transfer. The calculation is from a Heisenberg spin wave model 
with a small damping parameter (r=3 meV) (see text). 



from a standard vanadium sample was used to normalize 
the data so that S{Q^uj) is reported in units of mbarn 
sr~^ meV~^ per formula unit of CaFe2As2 throughout. 
We used the MSLICE program [U to visualize the data 
and to take one and two dimensional cuts through main 
crystallographic symmetry directions for subsequent data 
analysis with TOBYFIT [l5| . Symmetry equivalent cuts 
and slices were added to improve statistics. We used 
TOBYFIT to simulate the scattering cross-section for 
damped antiferromagnetic Heisenberg spin waves convo- 
luted with the instrumental resolution function for com- 
parison with the experiment. We note that orthorhombic 
twinning is present and accounted for in the simulations. 

Fig. [T] shows a projection of the Ei = 450 meV data 
along the iY-direction and near the (lOL) antiferromag- 
netic Bragg peaks as compared to the structure fac- 
tor calculated from a Heisenberg spin wave model with 
small damping (within instrumental precision) using ex- 
perimentally determined parameters (discussed below). 
The sizeable exchange interaction between planes leads 
to variations in the structure factor along L. The data 
indicate magnetic excitations beyond 150 meV and ex- 
tending perhaps as high as 200 meV. Excellent agree- 
ment between model and data was achieved though with 
considerable spin wave damping for energies beyond 100 
meV as will be discussed later. 

Fig. [2] shows constant energy slices taken in the HK- 
plane at specific energy transfers such that L = 2j-\-l [i.e. 
centered at Qafm = (1, 0, L = 2j + 1)]. The data clearly 
show rings of scattering as expected for constant energy 
slices through a conical dispersion surface emerging from 
QafM' Ring-like features are well-resolved up to ~ 100 
meV [see Fig. [S^c)]. Beyond 100 meV [Fig. E^d)], the 
data show smeared out and broadened features centered 
at Qafm- Fig. [3] emphasizes this observation by show- 
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FIG. 2: Constant energy slices through spin waves in CaFe2As2 at 
T = 10 K as observed on MAPS (left) (intensity shown in absolute 
units from background level) compared with calculations using the 
fitted exchange parameters in the text (right) with small damping 
(r=3 meV) (see text). 



ing several cuts along H at energies corresponding to 
L = 2 j + 1 . Well separated spin- wave peaks are observed 
between 60 and 118 meV. Above 100 meV, splitting is 
however not clearly observed due to weakening and/or 
broadening of the scattering. 

The data cuts, such as those in Fig. [3l were fit to Gaus- 
sian lineshapes in order to determine the relative position 
of the peaks with respect to the BZ center. The peak po- 
sitions were used to construct the spin wave dispersion 
relation along major symmetry directions (Fig. [H) These 
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FIG. 3: Intensity of observed spin wave excitations in CaFe2As2 
along [ifOG] direction at T = 10 K, compared with the resolution 
convoluted spin-wave calculations including damping shown in Fig. 
5. The calculations with energy independent damping (F held fixed 
at 3 meV) are also shown for comparison. 



data are compared to the dispersion relation in the lin- 
ear approximation for a Heisenberg model with nearest- 
neighbor (NN) and next-nearest-neighbor (NNN) inter- 
actions. The associated dispersion relations are given by 



E{ci) = JAl-Bl, 



(1) 



where 



= 2(5')[Ji5(cos(7rA;) 

5q = 2{S)[JlaCOs{^h)- 

+Jic cos(7r/)]. 



- 1) + Jla + Jlc + 2 J2 

2 J2 cos(7r/i) cos(7r/c) 



D] 



(2) 



Here the NN Heisenberg exchange constants are denoted 
Jiai Jihi and Jlc and the in-plane NNN exchange con- 
stant J2. {S) is the ordered spin, q = (hkl) = Q — Qafm 
is the reduced wave vector, with Q the wave vector 
transfer to the sample, and is a single ion uniaxial 
anisotropy constant. The corresponding equations for 
the spectral weight of the spin waves used in the simula- 
tion together with Eq. [T]are given explicitly in 

Heavy damping at high energies renders the Brillouin 
zone boundary spin wave energies ill defined. As a re- 
sult, only three independent energy scales can be ex- 
tracted from the data and these are Ji = (Jia + Ji5)/2, 
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FIG. 4: The dispersion of magnetic excitations in CaFe2As2: the- 
ory (lines) and experiment (solid circles). The dashed black lines 
are from Eq. fusing exchange parameters derived from band struc- 
ture calculations (see Ref. [lO]). The solid red lines represent the 
best global fit of Eq. [T] to the data. Variations of the fits with 
changing BZ boundary energy are shown for Ebz=^0 and 150 
meV. Also shown with open circles are the location of maxima in 
Im{x (q, <-t;)} from ab-initio band structure results. 



J2 = ^2 + ^^-^1/2, and Jlc- Here SJi = ( Jia — -^i6)/2, is re- 
lated to the putative zone boundary energy, which when 
Jlc << Jla and Jlc « J2 is given by £^(q = (010)) = 

4:{S)\j25Ji{2J2 — Ji). Fits with full resolution convolu- 
tion to the cuts in Figs. 1-4 yield {S)Ji = 21.9(9) meV, 
{S)J2 = 34.8(4) meV, and {S)Jic = 4.5(1) meV. We 
maintained the zone center gap at A = 7 meV in accor- 
dance with our higher resolution low energy data[lo| and 
this corresponds to {S)D ^ (AV8(^)( Ji + 2 J2 + Jic)) = 
0.063 meV. We obtain, g{S)iiB = 0.8(l)/iB from the in- 
tensity and this is consistent with the observed ordered 
moment. If 50 meV< £^(q = (010)) < 150 meV these re- 
sults indicate that 2 meV< {S)SJi < 15 meV such that 
24 meV< {S)Jia < 37 meV, 7 meV< {S)Jib < 20 meV, 
and 28 meV< {S)J2 < 34 meV. The fitted values are 
consistent with the requirement that {S){Jia + Jib) = 
2(6') Ji = 44 meV< 4{S)J2 = 124 meV±12 meV for the 
'stripe' AFM structure. The results also confirm that 
{S)Jc = 4.5(1) meV is substantial and the magnetic sys- 
tem is three-dimensional (our previous measurements [lO| 
could only establish a lower bound for {S)Jc). Fig. [H 
summarizes the fit results. 

The fits described above do not require substantial 
damping at low energies and the calculated intensities 
agree reasonably well there with no adjustable scale fac- 
tor (see Fig. [3|). This indicates that the full static 
moment participates in long wavelength collective (spin 
wave) excitations. To account for the data beyond 100 
meV it is however, necessary to introduce considerable 
spin wave damping. For simplicity the spin wave relax- 
ation rate, T{E)^ was assumed to depend only on en- 
ergy and was treated as piecewise constant in 10 energy 
bands. The simulations shown in Fig. 3 include energy 
dependent damping required to fit the data and Fig. 5 
shows the corresponding r(£^). This analysis provides ev- 
idence for Landau damping in the long wavelength limit 
and strong intrinsic broadening beyond 100 meV. Both 
findings are consistent with expectations for an itiner- 
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FIG. 5: Spectral half width at half maximum derived from band- 
structure calculated peaks in Im{x (q, a;)} (open circles) compared 
to the spin- wave relaxation rate, T{E) required to fit the neutron 
scattering data (solid circles). 

ant antiferromagnet. The situation is remarkably similar 
to underdoped YBa2Cu3 06.35 {Tc = 18 K) where there 
is an onset of damping and reduced spectral weight as 
compared to conventional spin wave theory beyond the 
150-200 meV pseudo-gap [17]. 

To better understand the origin of this broadening, 
theoretical calculations of the dynamic spin susceptibil- 
ity X (q, cj) were performed using linear response density 
functional approach [l^]. The calculated magnetic mo- 
ment on the Fe atom is 0.95/iB- Fig- [H shows a remark- 
able agreement between the calculated spin wave spec- 
trum, as derived from the peak position of Im{x (q, cj)} 
along the iY-direction, the raw experimental peak posi- 
tions (solid circles), and the fitted dispersion relation. 
The half- width at half maximum of the calculated peaks, 
are compared to the spin-wave relaxation rate extracted 
from the scattering data in Fig. [5l Strong line broad- 
ening is apparent beyond a; ~ 100 meV (or reduced 
g(lOO) ~ 0.2) and persists up to a maximum spin wave 
energy of 200 meV where T/co approaches 25%. The 
calculations indicate the presence of incoherent spin fluc- 
tuations across most of the BZ consistent with the ex- 
perimental data. Numerous band structure calculations 
indicate that the Fermi level sits in a region of relatively 
low density of states (pseudo-gap) approximately 0.1-0.2 
eV wide. The pseudo-gap may serve as a source of sta- 
bility for local moment collective excitations occurring 
at small q < 0.2 where a Heisenberg model description 
is justified. At higher energies, ab-initio calculations are 
best described as gradual evolution of the spin dynamics 
into the Stoner continuum. 

The present results favor an itinerant model for mag- 
netic excitations in the parent CaFe2As2 compound, as 
opposed to a local moment Heisenberg model. At small 

an effective NNN Heisenberg model is always an appro- 
priate description for either a local moment or itinerant 



system. The real applicability of such a model should be 
tested at larger q {q > 0.2). At these q both experiment 
and theory show the appearance of substantial damp- 
ing. Thus, the parent compounds already manifest strong 
coupling of magnetism and charge carriers in the form of 
an evolution of Landau damping to a full particle-hole 
continuum. This general remark should also apply to the 
doped superconducting compositions where AFM order- 
ing is suppressed but essential AFM correlations persist 
up to high energies [3]. In these compositions, it appears 
that strong Landau damping of the magnetic excitations 
occurs at much lower energies (~ 10 meV) and the SC 
gap suppresses damping. 

In summary, magnetic excitations in antiferromagnet- 
ically ordered CaFe2As2 have been observed up to 200 
meV. The features of the highly dispersive spin waves 
are consistent with local moment magnetism for energies 
below ~ 100 meV and can be represented by a NNN 
Heisenberg model there. For energies above ~ 100 meV, 
the excitations are strongly damped indicating presence 
of particle-hole excitations. CaFe2As2 is thus best de- 
scribed as an itinerant three-dimensional antiferromag- 
net. 

Work at the Ames Laboratory and the Johns Hop- 
kins University was supported by the Department of 
Energy, Basic Energy Sciences under Contract No. 
DE-AC02-07CH11358 and DE-FG02-08ER46544 respec- 
tively. Technical assistance of ISIS staff is gratefully ac- 
knowledged. 



[1] I. Mazin and M. Johannes, Nature Physics 1160, (2008). 

[2] Y. Qiu et a/., Phys. Rev. B 78, 052508 (2008). 

[3] K. Matan et a/., arXiv (2008), 0810.4790vl. 

[4] A. D. Christiansen et al , Nature 456, 930 (2008). 

[5] M. D. Lumsden et al, arXiv (2008), 0811.4755vl. 

[6] S. Chi et al, arXiv (2008), 0812.1354. 

[7] I. Eremin et al, Phys. Rev. Lett. 101, 187001 (2008). 

[8] A. V. Chubukov and L. P. Gor'kov, Phys. Rev. Lett. 101, 

147004 (2008). 
[9] J. Zhao et al, Phys. Rev. Lett. 101, 167203 (2008). 
[10] R. J. McQueeney et al, Phys. Rev. Lett. 101, 227205 

(2008). 

[11] R. Ewings et al, Phys. Rev. B 78, 220501 (R) (2008). 

[12] N. Ni et al, Phys. Rev. B 78, 014523 (2008). 

[13] A. I. Goldman et al, Phys. Rev. B 78, 100506 (R) (2008). 

[14] R. Coldea, MSLICE, A Data Analysis Programme for 
TOF Neutron Spectrometers (2004). 

[15] T.G. Perring, TOBYFIT version 2.0, Least-Squares Fit- 
ting to single crystal data (2004). 

[16] N. E. Zein and V. P. Antropov, J. Appl. Phys. 89, 7314 
(2001), Phys. Rev. Lett. 89, 126402 (2002) 

[17] C. Stock et al, Phys. Rev. B 75, 172510 (2007). 



